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Abstract. Possible consequences of ferromagnetic transi- 
tion in dense matter suggested recently by Kutschera and 
Wojcik, for the magnetic properties of neutron stars, are 
studied. Specific model of dense matter, in which a small 
admixture of protons is completely polarized due to their 
interaction with neutrons, is considered. Magnetic field 
of neutron stars with a ferromagnetic core are calculated 
within the framework of general relativity. Two types of 
boundary conditions at the ferromagnetic core edge are 
considered, corresponding to normal and superconduct- 
ing liquid envelope, respectively. Numerical results for the 
neutron star magnetic dipole moment are confronted with 
pulsar timing. To be consistent with observations, ferro- 
magnetic cores surrounded by a non-superconducting en- 
velope, should consist of weakly ordered ferromagnetic do- 
mains. If domains are highly ordered, ferromagnetic core 
should be screened by a superconducting envelope. 
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1. Introduction 

Observations of radio and X - ray pulsars indicate, that 
the surface magnetic field of neutron stars can be as high 
as 10^^ — 10^^ G. The strength of the internal magnetic 
field of neutron stars is unknown; it could be much higher 
than that of the surface one. 

Standard present day models of magnetic field iden- 
tify its source with long-lived electric currents fiowing in 
highly conductive neutron star matter. However, a com- 
plete scheme, based on such a model, which would explain 
both the origin and the evolution of neutron star magnetic 
field, is still lacking (recent review on the present status 
of this topic can be found in Bhattacharya & Srinivasan 
1995). 
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A very different type of model of neutron star mag- 
netic field is based on the hypothesis of existence of a 
ferromagnetic core in the liquid interior of neutron star. 
First models of this type were proposed just after the dis- 
covery of pulsars (Brownell & Callaway 1969, Rice 1969, 
Silverstein 1969). Neutron star matter, approximated by a 
pure neutron matter, was assumed to undergo a transition 
to a ferromagnetic state above some critical density. Fer- 
romagnetic neutron matter, carrying high magnetization 
density, was there a permanent source of a superstrong 
magnetic field. This field, which would permeate stellar 
interior, was expected to extend outside the stellar sur- 
face as an observable (mostly dipole) magnetic field of a 
pulsar. The qualitative argument in favor of a ferromag- 
netic transition in sufficiently dense neutron matter was 
based on Pauli exclusion principle, combined with the re- 
pulsive character of the short range n — n interaction. At 
high density, the spin singlet {S — 0) n — n interaction 
becomes strongly repulsive, while the S = 1 one, due to 
Pauli principle, avoids the most repulsive, shortest range 
contribution from the I = state (only odd / states are al- 
lowed in the S' = 1 channel). So, although in the ferromag- 
netic transition the kinetic energy of the system increases, 
at sufficiently high density this was believed to be more 
than balanced by the removal of the repulsive interactions 
in the S = states (because of complete spin polariza- 
tion, all neutron pairs in a ferromagnetic phase are in the 
spin triplet state). The argument seemed to be particu- 
larly convincing for the schematic n — n interaction of the 
infinite hard-core type (hard spheres gas). Some early cal- 
culations, based on the schematic models of the n—n inter- 
action, seemed to show the existence of spontaneous spin 
polarization in dense neutron matter above some critical 
density (Brownell & Callaway 1969, Rice 1969, Silverstein 
1969, Ostgaard 1970; see, however, Clark & Chao 1969, 
Clark 1969). However, further calculations, based on more 
realistic n — n interactions and/or more precise methods of 
solution of the many body problem, ruled out possibility of 
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ferromagnetism of dense neutron matter (Pandharipande 

ct al. 1972, Haensel 1975, Modarres & Irvine 1979). Such 
a possibility was ruled out even in the most favorable case 
of the hard-sphere model of neutron matter, where the 
contribution from the I = 1 states was shown to be very 
important at superhigh densities (D§,browski et al. 1978, 
1979). 

It should be mentioned, that in contrast to the realistic 
calculations described above, the calculations done within 
the Hartree-Fock approximation using schematic effec- 
tive nucleon Hamiltionian (or Lagrangian) often showed 
a ferromagnetic transition in neutron matter at a few 
times nuclear matter density. Vidaurre et al. (1984) used 
non-relativistic effective n — n interaction of the Skyrme 
type, and found a ferromagnetic transition at about 1.5 
times nuclear density. However, their result might be just 
a consequence of the specific density dependence of the 
spin dependent terms in the Skyrme interaction, which in 
most cases leads to unphysical collapse (no lower bound 
for the energy density) of dense polarized neutron mat- 
ter (Kutschera & Wojcik 1994). The spin stability of 
dense neutron matter was also studied within the relativis- 
tic Dirac-Hartree-Fock approximation, with an effective 
nucleon-meson Lagrangian, by Marcos et al (1991). These 
authors find a ferromagnetic transition at several times 
nuclear density, and point out that the presence of isovec- 
tor mesons in the mean-field Lagrangian is crucial for the 
possibility of ferromagnetic transition. Let us notice, that 
the importance of vector meson-coupling for the ocurrence 
of ferromagnetism was pointed out some 20 years earlier 
(Kalman & Lai 1971). In general, it seems that the fer- 
romagnetic transition in the mean-field theory models of 
dense neutron matter is a consequence of a very specific 
density dependence of the effective interactions character- 
istic of these models. 

In all studies reviewed above, neutron star matter was 
approximated by pure neutron matter. The importance 
of the presence of a small admixture of protons in neu- 
tron star matter for the spin stability was pointed out by 
Kutschera and Wojcik (1989). They assumed that neu- 
trons and protons in neutron star matter form normal 
Fermi liquids. The ferromagnetic state, corresponding to 
completely polarized protons and weakly polarized neu- 
trons, was shown to be energetically preferred over the 
non-polarized one in two cases: for sufficiently strong 
spin-spin component of the neutron-proton quasiparticle 
interaction and/or for localized protons. Such conditions 
were shown to be satisfied for numerous models of the 
neutron star matter with low proton fraction. With their 
models of ferromagnetic cores in neutron stars, Kutschera 
and Wojcik were able to make specific predictions, con- 
cerning dependence of the strength of the surface mag- 
netic field on the neutron star mass (Kutschera & Wojcik 
1992). The calculations performed for strongly asymmet- 
ric nuclear matter within the relativistic Dirac-Hartree- 
Fock approximation, using the mean field nucleon-meson 



Lagrangian, confirmed that the presence of an admixture 
of protons favors the ferromagnetic instability of dense 
matter (Bernardos et al. 1995). 

At the densities exceeding the critical density for the 
onset of ferromagnetism, the ground state of neutron star 
matter corresponds to a non-zero spin polarization. How- 
ever, the direction of the local spin polarization could vary 
within the ferromagnetic core, which might have a do- 
main structure. The macroscopic magnetization, relevant 
for the stellar magnetic field, may be expected to result 
from some volume averaging, and should be consistent 
with hydrostatic equilibrium of the star. We present in 
this paper a self-consistent calculation of the magnetic 
field of neutron star, produced by a liquid, ferromagnetic 
core. Our calculations, which are done within the frame- 
work of general relativity, take into account constraints 
resulting from the requirement of the hydrostatic equilib- 
rium, as well as from realistic boundary conditions. 

The paper is organized as follows. In Section 2 we 
present the microscopic model of nucleon ferromagnetism. 
A model of ferromagnetic neutron star core is described in 
Section 3. Physical conditions within the neutron star in- 
terior, relevant for the properties of ferromagnetic neutron 
stars, are discussed in Section 4. The calculation of mag- 
netic field of neutron star with a ferromagnetic core, per- 
formed in general relativity, is described in Section 5. Nu- 
merical results, obtained using specific assumptions con- 
cerning physical conditions in neutron star interiors, are 
presented in Section 6. In Section 7 we confront our re- 
sults with observations of radio pulsars. Our conclusions 
are presented in Section 8. 

2. Ferromagnetism of dense matter 

The possibility of ferromagnetism of dense nucleon mat- 
ter results from the existence of the spin-spin component 
of the quasiparticle nucleon-nucleon interaction. We con- 
sider the simplest model of neutron star matter, consist- 
ing of neutrons, protons and electrons. At given baryon 
density, rib, the composition of matter can be character- 
ized by the proton fraction x, such that proton density 
Up = xrib = rie, while neutron density rin = (1 — x)nh. 
The degree of spin polarization of nucleons can be mea- 
sured by the parameters 



where rini is the number density of neutron quasiparticles 
with -1-1/2 spin projection of the spin on the spin quan- 
tization axis, etc. In the normal (i.e., non-ferromagnetic) 
state, the energy density reaches its minimum at s„ = 
Sp = 0. 

Our discussion of the properties of the spin-polarized 
neutron star matter will follow the paper of Kutschera 
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and Wojcik (1989). We assume that nucleons arc normal 
(possible effects of nucleon superfluidity will be discussed 
in Sections 4, 6). We neglect the thermal effects, and cal- 
culate all thermodynamic quantities at the T = approx- 
imation. The change in the nucleon contribution to the 
energy density, implied by a small spin polarization of nu- 
cleons, can be calculated using the methods of the 
theory of Fermi liquids. 
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where Gq", and g^^ are Fermi liquid parameters, de- 
scribing the spin-spin terms in the qusiparticle nucleon- 
nucleon interaction in nucleon matter, and N^, Np are 
the densities of states at the corresponding Fermi sur- 
faces for neutrons and protons, respectively. While first 
and and second terms on the r.h.s. of Eq.(2) are posi- 
tive (Gq", Gq^ > — 1), the last term can be always made 
negative, irrespectively of the sign of the ^q'' parameter. 
Notice, that the third term, which results from the spin- 
spin interaction between neutron and proton quasiparti- 
cles, vanishes for a non-interacting system. 

If (5£'^i < 0, the system is unstable with respect 
to spin polarization. For small proton fraction {x <C 1), 
the approximate expression for SE^^^ reads (Kutschera & 
Wojcik 1989, hereafter referred to as KW89) 
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Spontaneous spin polarization takes place for sufficiently 
strong Qq^. The approximate condition for nucleon ferro- 
magnetism reads thus 



(97) > 



1 



Gnn 




NnNp 



(4) 



If condition (|j) is satisfied, the ground state of matter 
contains then completely polarized protons (sp = 1), and 
weakly polarized neutrons. 



l + Gg" 



(5) 



The values of the Fermi liquid parameters in strongly 
asymmetric nuclear matter, as well as their density de- 
pendence, are very uncertain. It seems to be reasonable 
to say, that their present knowledge does not exclude the 
situation, in which Equation (^) is satisfied above some 
critical density n*. 

For sufficiently low proton fraction (such that protons 
can be considered as impurities in neutron matter), the 
ground state of nucleon matter could correspond to a 
non-uniform distribution of neutrons, with proton local- 
ized (bound) in the neutron density minima (Kutschera & 



Wojcik 1990). In the case of localized protons, the whole 
second term on the r.h.s. of Eq.(2) vanishes, and nucleon 
matter is always unstable with respect to full polarization 
of protons (KW89). Such a case corresponds most proba- 
bly to the solid ferromagnetic core (Kutschera & Wojcik 
1995). 

It should be stressed, that the results described above 
can be directly used in the neutron star calculations only 
under a rather unrealistic assumption of perfect space ho- 
mogeneity of a macroscopic element of matter. Only in 
such an idealized case (which was actually assumed in 
Kutschera & Wojcik 1992) the 'microscopic' and 'macro- 
scopic' magnetizations will coincide. It seems to be quite 
realistic to assume, that the ferromagnetic phase of dense 
matter has some 'domain structure'. In such a case, micro- 
scopic results would be valid only within a single domain. 
The macroscopic values, relevant for the stellar structure, 
would then correspond to an average over a sufficiently 
large number of domains. To be more precise, the macro- 
scopic value of a quantity is then defined as an average 
over volume element, which contains a sufficient number 
of domains, so that the calculated value does not depend 
on the specific averaging procedure. Of course, the ele- 
ment should be still small on the stellar scale, so that the 
macroscopic value has a local character. The macroscopic 
quantities would have to be consistent with our assump- 
tion about the axial symmetry of the star. 

It is well known, that ferromagnetism is an example 
of a second-order phase transition, and can be treated 
within the general scheme of the theory of second-order 
phase transitions, formulated by Landau (see, e.g.. Sec- 
tion 39 of Landau et al. 1984). However, Landau theory 
of second-order phase transitions asumes expansion of E 
(or other suitable thermodynamic potential) up to fourth 
order in polarization parameter. Approximate expression 
for SEp^^, Eq. (2), is truncated at terms quadratic in po- 
larization parameters. Such an approximation is charac- 
teristic of the Fermi liquid theory, in which only terms 
quadratic in deviations from the normal, reference state, 
are conserved. Therefore, the model proposed in KW89, 
if taken literally, cannot predict the threshold behavior 
Sp, Sn c>c (n — n*)^/^ for n — > n» -1- 0, characteristic of a 
general case of spontaneous polarization in ferromagnetic 
phase. 

The case of ferromagnetism with localized protons de- 
serves a special comment. The threshold density for proton 
localization, nioc, corresponds to a phase transition from 
a liquid mixture of protons, neutron and electrons to a 
proton crystal immersed in a neutron and electron liquid. 
Such a liquid-solid phase transitions is expected to be the 
first-order one (i.e., accompanied by the density and order 
parameter discontinuity). However, for normal neutrons, 
localization implies complete proton polarization, so that 
n* — n\oc- In such a case, even for rib — > + protons 
remain fully polarized. 
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3. Ferromagnetic neutron star cores 

For simplicity, we assume that magnetization of the fer- 
romagnetic core represents the only source of stellar mag- 
netic field (i.e., we neglect contribution to magnetic field 
resulting from electric currents). 

We consider a simple model for the magnetization den- 
sity of polarized matter, m, within the ferromagnetic core 
(rib > n,), as a function of baryon density, rib- In gen- 
eral, the spin contribution to the z-component of m in the 
ferromagnetic phase is given by 

(^"^)spin = npSpUp + n„s„^„ , (6) 

where ^p, fin are intrinsic magnetic moments of protons 
and neutrons, respectively, and we assumed that z-axis 
is the spin quantization axis. We modeled (™z)spin by a 
simple expression 

im.Kin-^^f", (7) 

where m™"^'' corresponds to perfectly ordered ferromag- 
netic state, considered by Kutschera & Wojcik (1992, here- 
after referred to as KW92), and dimensionless factor a < 1 
describes the degree of macroscopic polarization (degree 
of ordering of ferromagnetic domains). Our formula for 
m, Eq.(|^), takes thus into account a possibility of partial 
macroscopic ordering of nucleon spins (partial ordering of 
ferromagnetic domains), including the case of vanishing 
macroscopic magnetization within the ferromagnetic core 
(« = 0). 

Macroscopic magnetization density m is the source of 
magnetic field within ferromagnetic neutron star cores. On 
the other hand, the presence of non-vanishing m modi- 
fies the energy density, due to coupling to magnetic field. 
Consistent derivation of this coupling from microscopic 
theory shows, that the additional term to be included is 
— nispin • B (de Groot & Suttorp 1972, Chapter X; Carter 
1982). If we take the z-axis along B, then minimization 
of energy requires nizBz > 0. The consequences of this 
condition will be discussed in Section 6. 

4. Physical conditions in neutron star interior 

In what follows, we consider the interior of neutron star, 
of density rib > ?T.h = 0.1 fm~^. With our assumptions, 
this interior consists of a ferromagnetic core of the density 
rib > n^, and a non-ferromagnetic liquid envelope of the 
density nh < rib < n*. 

Neutron stars are believed to be born as very hot 
objects, with initial internal temperatures above 10^^ K. 
Young neutron star cools via neutrino emission. Initial 
temperature of neutron star core can be expected to be 
larger than, or of the order of, the Curie temperature for 
the ferromagnetic transition, Tferro ~ few times 10^° K 
(Haensel 1995, unpublished). Another important criti- 
cal temperature is that for nucleon superfluidity, Tsup ^ 



10^ K. It is thus reasonable to expect, that the nucleon fer- 
romagnetic transition occurs before superfluidity sets in, 
i.e., that ferromagnetism takes place in normal Fermi liq- 
uids. Because Tfeno > Tsup, completely polarized proton 
component cannot undergo singlet pairing at Tgup- Also, a 
fraction of polarized neutrons remains in the normal state, 
due to their coupling to completely polarized protons. In 
such a way, protons, and a small fraction (a few percent) 
of neutrons may be expected to be locked in a polarized 
state, due to energy barrier against the depolarization of 
the ferromagnetic component. 

It should be stressed, that in view of the uncertain- 
ties in the many body theory of dense matter at supranu- 
clear density, even the very existence of proton super- 
conductivity cannot be considered as absolutely certain. 
In view of this, we may contemplate several possibili- 
ties. In the first case, considered by KW92, ferromagnetic 
core is surrounded by a liquid envelope of normal mat- 
ter with nil < ^b < '^*, containing normal (i.e., non- 
superconducting) protons. In the second case, liquid en- 
velope is assumed to contain superconducting protons. 
The superconductivity of the envelope could strongly in- 
fluence magnetic properties of ferromagnetic neutron star 
(see Section 5, 6). 

5. Calculation of magnetic field 

The calculation of the stellar magnetic field has been done 
neglecting sources other than magnetization of matter 
within the ferromagnetic core. We took into account space 
curvature, implied by the neutron star gravity. 

The first stage consisted in the calculation of the con- 
figuration of hydrostatic equilibrium for a non-rotating 
neutron star with a ferromagnetic core. In this calcula- 
tion, we neglected the infiuence of the presence of mag- 
netic field on the equation of state of the matter; this 
approximation will be justified quantitatively in Section 

6. The magnetic field vectors H^, were calculated by 
solving the equations of magnetostatics in curved space. 
These equations are deduced from the Maxwell equations 
in the curved space in the static case (see, e.g., Eqs. (7.6) 
- (7.7) of Carter (1980), or Thorne & Macdonald (1982)), 

ViS' = , (8) 

e'j'^Vj(7Vi?'^) = , (9) 

where Vi denotes the covariant derivative on the t — const 
hyperspaces, N is the lapse function, which in the static 
case is given by iV = \/— 5ooj and e'j'^ is the antisymmetric 
Levi-Civita tensor associated with the 3-metric 7ij on the 
t = const hypersurfaces. Equation (||) will be satisfied, if 
and only if, is the gradient of some scalar function 

^, so that 

- i-ViVi/ . (10) 
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Inserting ( |10[ ) into (g) we get, using = i?' + 47rm' 
1 



A^f = -47rViTO' 



N 



(11) 



where A is the Laplacian operator in the space curved by 
neutron star gravity, given by 



1 



(12) 



where 7 = det(7ij). In view of the smallness of the defor- 
mation of the star, implied by the presence of magnetic 
field, we can consider the non-trivial elements of the met- 
ric 7ij as depending on the radial coordinates only. 



711 = A'{r) , 722 =rM'(r) 



733 



r sm 
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Within this approximation, equation (12) can be rewritten 
as 



A^- 



1 d 
r^A^ dr 
1 f 92 



or 
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In what follows, we will make a simplifying assumption 
that m(r) is parallel to the symmetry axis of the star. 
The non-vanishing normal components of magnetization 
in the local coordinate system are then: = to cos 6*, 
mg = —msmO (where to = |m| depends on the r coordi- 
nate only), so that toi = Anir, m} = rrir/A, — Arrag, 
TO^ = mg/{Ar). The divergence VjTO' reads, in our ap- 
proximation, 

„ i cos 6 d 

ViTO' = 



{A^ 



A^ dr 

The Eq. ( pr] ) can thus be rewritten as 
1 d f ^^d \ A f ^ 



(15) 
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dr 



-ATTCOsd— (A 
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00^ 
A dN a* 
N dr dr 
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Let us remind, that within our approximations A, N , and 
TO are functions of r only. Putting 'J — <I>(r)cos0 we 
rewrite Eq. (16) in a more suitable form: 



1 
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dr 



r^A^<^ 
dr 
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-4^— (A^ 
dr ^ 

A dN d$ 

N dr dr 



(17) 



As a consequence of general properties of the second- 
order phase transitions, the derivative of to with respect 
to r is singular at the edge of ferromagnetic core, r — r^. 
This is implied by the behavior to — > const. (nt, — r!*)^/^ 



for rih — > -f 0. To avoid numerical problems stemming 
from this singularity, we rewrite Eq. (17) as 



d 
dr 



Ar^^(rH) 
dr ^ ' 



-47r-^ (A^m) 
dr ^ ' 

2 d , A diV d$ 

-$ — A H 

r dr dr dr 



(18) 



which implies 



Ad_ 

j,2 



d ,^ „dAd$ 
2^-^(0+A(r')-- 



dr' 



dr' 



(19) 



so that 



$(r) = -% I r''^A^{r')m{r')dr' 

T Jo 



dr' / X 



1 

Jo A{r')"' Jo 
A<i>(r")^A(r")+A(r")^^ 



^dA^ d$ 
dr" dr' 



dr" . (20) 



Let us consider first the case of normal envelope of the 
ferromagnetic core. Then, Equation (20) can be solved 
by iterations in all the space, < r < 00, with regu- 
larity condition ^(O) = 0, and with boundary condition 
^(cxd) = 0. Solutions were obtained using pseudospectral 
method, with three spatial grids, in each of the three re- 
gions 0<r<r^, r» < r < R and R < r < 00, respec- 
tively. The method used for compactification of the space 
is the same as that used for the calculation of the met- 
ric gij (see Bonazzola et al. 1993). Once $ is known, the 
vectors and are calculated from 



Hr{r, 6) 
Heir, 6) 
Br{r,e) 
Bg{r, 9) 



1 /. X 

- — <i> (r) cos( 
AN ^ ' 



ANr 



$(r) sinfl 



= — — <i>'(r) COS0 -I- 47rTOr , 

1 

<i>(r) sin 9 + Anmg 



ANr 



(21) 



Magnetic permeability of normal neutron star matter is 
close to unity; for simplicity, we replaced it by one. 

In the second case, we assume that the ferromagnetic 
core is surrounded by liquid superconducting envelope 
(with superconducting protons). Here, we can deal with 
two distinct situations. If protons form a type-I supercon- 
ductor, then B = for nh < nb < (Meissner effect). 
This would correspond to a perfect screening of the fer- 
romagnetic core. It seems more likely, however, that the 
magnetic field nucleated within the superconducting enve- 
lope as an array of fluxoids, a situation characteristic for 
a type-II superconductor (Baym et al. 1969, Sauls 1989). 
Such a superconducting shell would behave as a diamag- 
netic envelope of ferromagnetic core. 
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6. Magnetic field under various physical conditions 

Our models of neutron star were calculated assuming the 
BJI model of the equation of state (EOS) of neutron star 
matter (Malone et al. 1975). With this EOS, we obtain 
neutron star models which are quite similar (for M 1.4— 
1.8 Mq) to those calculated using more recent EOS, such 
as the UV14+TNI or FPR models considered by KW92. 
For neutron star mass M — 1.41 M©, we get the radius 
R = 12.11 km, and central density nccntr — 0.58 fm~'^. 

In our calculations of magnetic field, we assumed con- 
stant value of the degree of macroscopic polarization, 
a, and a constant, small value of the proton fraction, 
X = 0.05. In view of linearity of Eq. (17), the values of B, 
H for a specific value of a can be obtained from the a = 1 
ones using B — aB{a = 1), H = aH(Q! = 1). For the 
model of KW92, the microscopic magnetization m™"^'' has 
a specific dependence on the nucleon density, rib = rin+Up. 
We fitted the density dependence of m™^" by an analytic 
formula, which exhibited correct threshold behavior for 
nb ~ n*. We assumed — 0.3 fm~^. The specific density 
dependence of m™"^'' at intermediate densities results from 
the density dependence of Fermi liquid parameters, com- 
bined with opposite signs of fin and fip, as well as from the 
assumed dependence of the proton fraction on the nucleon 
density. The quantity m™"^'', calculated in KW92 using Eq. 
(^, with Sp = 1 and Sn given by Eq. (^, changes sign at 
some rii-. However, the condition mB > puts a definite 
condition on the sign of ruz on the stellar symmetry axis, 
which reads m^Bz > 0. In view of this, vanishing of at 
some point (with ^ 0) should not be accompanied by 
the change of sign of in the vicinity of this point. In 
view of this, passing through n,- does not change the sign 
of rUz, but requires the reversal the signs of Sn and Sp. Let 
us notice, that the nuclear contribution to the polarization 
energy, (5i?^j, is invariant with respect to the change of the 
orientation of spins, Sn — > — Sn, Sp — > — Sp, while such 
an inversion of spins changes the sign of m™"^''. Therefore, 
imposing the constraint m™"^''i?2 > in the case of a = 1, 
when iriz = m'^^", has no effect on the nuclear part of 
the energy density. However, it will have important conse- 
quences for the magnetic properties of the ferromagnetic 
core. 

Our basic calculations have been performed, assuming 
a — 1 and M = 1.41 Mq. We considered the case of a 
liquid envelope with normal protons. The radial distribu- 
tion of niz is displayed in Fig. 1. In Fig. 2 we show the the 
values of Hz and Bz, respectively, on the symmetry axis, 
as functions of z. 

Let us notice, that relativistic effects in the ferromag- 
netic core are quite large. From equations (19) and (20), 
we obtain relation between the values of Bz,Hz and the 
lapse function, N , at the center of the star. 



For our neutron star model we get Bz{0)/Hz{Q) = 
— 1.17, to be compared with ffat space-time value 
{Bzi'd)/ Hz{Q))f^^^ = —2 (M. Kutschera, private communi- 
cation) . The importance of relativistic effects results from 
the fact, that the lapse function at the star center is sig- 
nificantly lower than unity, N{Q) — 0.618. 
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Fig. 1. Magnetization of matter within the core versus distance 
from the star center, for the 1.41 Mq model. Perfect ordering 
of ferromagnetic domains ( q = 1 ) has been assumed. 



Let us notice, that the value of Bz at magnetic poles 
(z — zLR) is significantly smaller (~ four times) than its 
value at the ferromagnetic core edge. The dipole magnetic 
moment of the star is — 15.1 ( here Aisa is the 

magnetic dipole moment in the units of 10"^" G cm"^). At 
fixed M, the values of Bz{pole) and Ai^°"° for different 
values of the spin ordering parameter a can be calculated 
by multiplying by a the values obtained for a = 1. 

Both Bz{pole) and Ad^'^"" depend on the mass (cen- 
tral density) of the neutron star model. Obviously, we have 
"■contr < However, with our choice of 

= 0.3 fm~'^ and with our EOS, this implies that ferro- 
magnetic core is present in neutron stars with M > M^ = 
0.7 Mq - a condition likely to be satisfied by radio pul- 
sars. Both J\4^°™ and i3z(pole) increase with increasing 
stellar mass. We have A^so = 11.8, -Bz, 12 (pole) = 11.7 for 
M = 1.2 Mq, and both quantities icrease with increasing 
mass up to 7W30 = 108 and Bz,i2 — 174 for the maxi- 
mum allowable mass, M„iax = 1-85 Mq. Here, B12 is the 
magnetic field induction in the units of 10^^ G. 

It should be stressed, that in all cases the deformation 
of neutron star due to the presence of magnetic field is in- 
deed negligibly small. For our model, the maximum energy 
density of magnetic field created via the ferromagnetic 
magnetization is comparable to the gain in the internal nu- 
clear energy density resulting from the ferromagnetic tran- 
sition. A simple estimate, based on the formulae of Sect. 
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Fig. 2. The values of Hz (dotted line) and Bz (solid line) (in 
10^^ G) on the z-axis, as a function of the distance from the 
star center, for the 1.41 Mq star. Perfect ordering of ferromag- 
netic domains { a — 1 ) and normal liquid envelope have been 
assumed. 



2, gives 6E^^i « -3 10^\rn,/no)^/^ix/0m)^/^ erg/cm^, 
which is very small compared to the values of pressure 
near the neutron star core, P ^ 10^^ erg/cm'^. 



7. Confronting models with observations 

Results described in the preceding section should be con- 
fronted with existing information on the magnetic field 
of pulsars. When confronting ferromagnetic neutron star 
models with observations, one should be aware of the fact, 
that even if some pulsars do contain ferromagnetic cores, 
such a core would be most probably only one of the pos- 
sible sources of pulsar magnetic field, the other one being 
the long-living electric currents in the neutron star inte- 
rior. 

We will use a standard assumption, that the slow- 
ing down of pulsar rotation is due to the loss of its ro- 
tational kinetic energy, implied by the emission of the 
low frequency dipole radiation. We can then express the 
dipole magnetic moment a pulsar of period P and the time 
derivative of the period, P, as 



M 



PSR 



3cUPP 
87r2 sin^ /3 



(23) 



where f3 is the angle between the magnetic and rotation 
axis, and / is the moment of inertia of neutron star. As- 
suming canonical mass of radio pulsars, 1.4 Mq, which 
for our EOS corresponds to / = 1.45 x 10"*^ g cm^, we 
see that pulsar timing puts limits on the magnetic mo- 



ment of ferromagnetic core (expressed here in the units of 
1030 G cm3). 



Mf^'°\ sin I3\ < 14 



P 



10 ms 



P 



10 



-13 



(24) 



Available timing data from Taylor et al. (1993) yield a very 
large range of Alaol sin /3\: from 2.3 x 10"'^ for second most 
rapid millisecond pulsar PSR B1957-h20, to 0.83 for PSR 
B0154+61. Excluding the case of nearly perfect aligne- 
mcnt of magnetic and rotation axes, we find two situations 
consistent with observations. In the first case, ferromag- 
netic core is characterized by a low degree of macroscopic 
ordering, which corresponds to a = 10~^ — 10~^. In the 
second case, the ferromagnetic core is screened by the su- 
perconducting (type-I) envelope, and therefore constraint 
( p4[ ) is irrelevant there. The case, when protons in the liq- 
uid envelope form a type-II superconductor would be an 
intermediate one; observational constraint on A4^'^"°, Eq. 
(22), would then be weakened due to partial screening of 
ferromagnetic core. 



8. Discussion and conclusions 

In the ferromagnetic model of Kutschera and Wojcik 
(1992), consistency of the obtained values of surface mag- 
netic field with pulsar timing was obtained via cancella- 
tions of the contributions to the stellar dipole magnetic 
moment coming from the inner and outer regions of the 
ferromagnetic cores. In particular, the low values of B for 
millisecond pulsars were explained by assuming that the 
magnetic moment of their ferromagnetic cores was some 
three orders of magnitude smaller than that of ordinary 
pulsars. According to Kutschera and Wojcik (1992), such 
a situation was made possible due to the cancellations, re- 
sulting from the change of sign within the ferromagnetic 
core. This was possible only within a very narrow interval 
of neutron star masses (AM ^ 0.1 Mq). 

Consistent calculations of the magnetic field pro- 
duced by the ferromagnetic core of neutron star, per- 
formed in the present paper, show that the projection of 
the macroscopic magnetization onto the symmetry axis 
cannot change sign within the core. In view of this, con- 
sistency with pulsar timing can be achieved only under 
some specific conditions. If the liquid envelope surround- 
ing ferromagnetic core is non-superconducting, than the 
ferromagnetic phase has to be highly disordered, with the 
degree of macroscopic polarization varying from ~ 10"^ 
for young pulsars, to ^ 10^** for the millisecond pulsars. 
In the case, when the liquid envelope is a type-I super- 
conductor, the ferromagnetic core produces no external 
magnetic field, which could be then produced exclusively 
by the currents within the neutron star crust. 
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Recently, Kotlorz & Kutschera (1994) pointed out an 
interesting possibility of existence of ferromagnetic pion- 
condensed quark cores, within sufficiently massive neutron 
stars. Their results imply 30'''° (quark) ~ 10^. Such huge 
values of A^^''™ can be reconcilied with pulsar timing only 
if either pion-condensed quark core has very low degree of 
macroscopic spin ordering, a < 10~^, or if this core is 
screened by the superconducting nucleon envelope. 
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